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ABSTRACT 

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf 
spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations 
of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo 
(-^vir.Mw)- Despite Leo I's large Galacto-centric space velocity (200kms~^) and distance (261 kpc), 
we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark 
matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed 
position and velocity of Leo I strongly disfavor a low mass Milky Way: if we assume that Leo I is 
the least bound of the Milky Way's classical satellites, then we find that Mvir,MW > 10^^ Mq at 95% 
confidence for a variety of Bayesian priors on Mvir,MW- In lower mass halos, it is vanishingly rare 
to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found 
to be less bound than Leo I, this lower limit on Mvir.MW would increase by 30%. Imposing a mass 
weighted ACDM prior, we find a median Milky Way virial mass of Mvir.MW = 1-6 x 10^^ Mq, with a 
90% confidence interval of [1.0 — 2.4] x 10^^ Mq. We also confirm a strong correlation between subhalo 
infall time and orbital energy in the simulations and show that proper motions can aid significantly 
in interpreting the infall times and orbital histories of satellites. 

Subject headings: galaxies: individual (Leo I) - Galaxy: kinematics and dynamics - Galaxy: halo — 
Galaxy: fundamental parameters - methods: iV-body simulations 



1. INTRODUCTION 

Satellite galaxies orbiting the Milky Way (MW) have 
been a rich source of astrophysical information and cos- 
mological confusion. Individually, these objects probe 
galaxy formation at the lowest masses, luminosities, and 
metallicities. As a population, the Milky Way's dwarf 
satellites provide one of the few existing tests of cosmol- 
ogy on small scales. They also provide a sample of trac- 
ers that can be used to estimate the mass distribution of 
the Milky Way's dark matter halo. Their small number, 
anisotropic spatial distribution, and unknown tangential 
velocities are all potential sources of biases and misun- 
derstandings, however. 

In particular, Leo I has long been known to play on 
outsized role when using satellite galaxies to determine 
the mass of the Milky Way's dark matter halo (Mvir,Mw; 
Zaritsky et al. 1989; Fich & Tremaine 1991; Kochanek 
1996) because of its large distance and high radial ve- 
locity in the Galacto-centric frame. In the analysis of 
Watkins ct al. (2010), which incorporates line-of-sight 
velocity information for 25 Galactic satellites in deter- 
mining Afvir,MWj Leo I alone contributes approximately 
27% to the mass estimate (i.e., including Leo I leads to 
mass estimates that are 25-30% larger than those that 
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exclude Leo I). 

Incorporating Leo I into an estimate of Mvir,MW gener- 
ally requires the assumption that Leo I is indeed bound. 
The timing argument (Kahn & Woltjcr 1959), for exam- 
ple, assumes the MW and Leo I to be on an elliptical 
orbit, with Leo I having made a recent pericentric pass 
about the Milky Way. Since Leo I is, by assumption, 
bound to the Milky Way in such analyses, a large Milky 
Way mass is required to match Leo I's observed distance 
and radial velocity (e.g., Zaritsky et al. 1989; Li & White 
2008; Sohn ct al. 2012b). If Leo I is unbound, however, 
the timing argument as applied to Leo I will almost cer- 
tainly overestimate Mvir,MW- Leo I's extreme kinematic 
properties have therefore led to substantial debate about 
its status as a bound satellite of the Milky Way (Zaritsky 
ct al. 1989; Byrd et al. 1994; Sales et al. 2007; Sohn et al. 
2007; Mateo et al. 2008). 

Theoretical models of halo collapse provide a further 
basis for predicting satellite orbits in ACDM. In partic- 
ular, the self-similar secondary infall model (SSIM; Fill- 
more & Goldreich 1984; Bertschinger 1985) has guided 
much of the progress in understanding halo collapse and 
virialization. In the simplest version of the SSIM - one 
where only radial motions are considered - the matter 
flow around a slight overdensity in the early universe can 
be derived in a straightforward manner. Shells of matter 
around the overdensity initially expand with the Hubble 
flow. Those shells that are bound to the perturbation 
eventually reach a turn-around radius and re-collapse, 
with more bound shells collapsing earlier; the virial ra- 
dius of a halo is approximately equal to half of the present 
day turn-around radius in the SSIM. In this model, then, 
all accreted material must be bound; otherwise, it would 
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not collapse onto the overdensity. 

Although the physical picture of accretion in a hier- 
archical universe is somewhat different from the smooth 
SSIM, analyses of iV-body simulations have shown that 
phase space is often structured in a very similar manner 
to SSIM predictions (e.g., Ascasibar et al. 2007; Diemand 
& Kuhlen 2008). With numerical simulations, it is pos- 
sible to compute the orbits of halos falling into larger 
halos. Studies of large volume cosmological simulations 
generally find that unbound orbits are very rare, in accor- 
dance with the expectations of the SSIM: Wetzel (2011) 
found that less than 2% of orbits are unbound at all 
redshifts, while Benson (2005) found that only 0.3% of 
all orbits escape the host halo. These analyses typically 
adopt the point mass limit and therefore underestimate 
the potential energy of the orbit coming from the mass 
distribution exterior to the halos; accordingly, these un- 
bound fractions should be upper limits to the true un- 
bound fraction. While the works of Benson and Wetzel 
studied orbits of merging dark matter satellites across a 
broad spectrum of dark matter halo masses, Deason et al. 
(2011) and Di Cintio et al. (2012) focused on Milky Way- 
mass systems, again finding that the fraction of unbound 
satellites is quite low. In this paper, we show that the 
eight highest resolution simulations performed to date of 
Milky Way-size halos all have unbound fractions of less 
than 1 in 1000 (0.1%). 

In accordance with predictions of the SSIM, iV-body 
simulations show that halos falling into larger halos have 
a characteristic orbital energy that is roughly invariant 
with time: for host halos with M^h- ~ 10^^ Mq, this en- 
ergy corresponds to an orbital velocity of « 1.1 Kir at 
the virial radius, with a 1 ct scatter of 25% (Wetzel 2011; 
see also Benson 2005; Khochfar & Burkcrt 2006). Since 
halo virial masses grow with time, due to physical mass 
accretion or to a time-varying virial overdensity thresh- 
old (Diemand et al. 2007b; Diemer et al. 2012), subhalos 
accreted at early times will be somewhat more tightly 
bound than those accreted more recently. Once an in- 
falling halo has interacted substantially with the host, 
its orbital energy may decrease owing to processes such 
as dynamical friction; much less frequently, it may gain 
orbital energy through three-body interactions. Fluctua- 
tions in the central gravitational potential of the host due 
to, e.g., virialization or mergers may also affect subhalos' 
energies. Subhalos that have been part of their host for 
a long time will therefore exhibit a wide range of orbital 
energies, whereas recently accreted subhalos should oc- 
cupy a much narrower range in energy space (and should 
be less strongly bound). 

This general picture has been confirmed by Diemand 
& Kuhlen (2008), who show a strong relationship be- 
tween the number of orbits a subhalo has completed and 
its position in radial velocity phase space for the Via 
Lactea (VL) simulation (Diemand et al. 2007a), and by 
Rocha et al. (2012), who found a clear correlation be- 
tween orbital binding energy and infall time for subhalos 
in the Via Lactea II simulation (VL-II; Diemand ct al. 
2008). Similar relationships are generically expected in 
the SSIM, although this model does not account for en- 
ergy loss processes such as dynamical friction. Note that 
both in the SSIM and in cosmological simulations, satel- 
lites falling into dark matter halos for the first time do 
so on bound orbits: a first infall need not, and almost 



always does not, imply an unbound orbit. 

The three body interactions discussed in the previous 
paragraphs provide a means of boosting the velocities 
of some subhalos to above the local escape velocity, how- 
ever. This possibility has received attention in the litera- 
ture recently, often in the context of exploring "extreme" 
objects such as Leo I. Sales et al. (2007) emphasized the 
existence of dark matter satellites that have gained en- 
ergy through three body interactions with other, more 
massive satellites inside of their host halos. Ludlow et al. 
(2009) examined subhalo orbits around Milky Way ha- 
los as well and showed that many satellites travel to 
large apocentric distances, much larger than would be 
expected from their initial orbits, after gaining energy 
through the aforementioned three-body interactions dur- 
ing their first pericentric pass. Even after gaining en- 
ergy, these satellites typically remain formally bound to 
their dark matter host, although their apocenters can 
(significantly) exceed their original turn-around radius. 
These "ejected" subhalos always have very low mass at 
infall (0.1% of their host's mass or less; Ludlow et al. 
2009) and are usually found as the companion of a more 
massive infalling satellite (Sales ct al. 2007). In general, 
ACDM simulations predict that many satellites once well 
within a host's virial radius will spend significant time at 
larger radn (Mamon ct al. 2004; GiU et al. 2005; Wang 
ct al. 2009), a natural consequence of the highly eccentric 
nature of subhalo orbits (the median apocenter to peri- 
center distance ratio for subhalos is approximately 6:1; 
Ghigna et al. 1998; Diemand et al. 2007b) 

While the uncertainty surrounding Leo I's status as 
a bound satellite is a source of frustration, this does 
not diminish the importance of Leo I in establishing the 
dark matter distribution around the Milky Way. There 
are few kinematic tracers for D > 80 kpc (see, e.g., 
Gnedin et al. 2010; Deason et al. 2012), and deriving 
mass constraints from theses tracers (typically blue hor- 
izontal branch stars) requires assumptions about their 
spatial distribution and velocity anisotropics (Dchnen 
ct al. 2006; Deason et al. 2012). 

In a companion paper (Sohn ct al. 2012b; hereafter. 
Paper I), we presented the first determination of Leo 
I's transverse velocity, based on proper motion measure- 
ments using the unique astrometric capabilities of the 
Hubble Space Telescope (HST). This provides a direct 
probe of the mass distribution of the Milky Way at 260 
kpc, where the potential is dominated by the dark matter 
halo, and therefore a new window on Mvir,MW- In this 
paper, we combine the measurement of the proper mo- 
tion of Leo I from Paper I with state-of-the-art numerical 
simulations of the formation of Milky Way-size halos in 
ACDM in order to constrain the mass of the Galaxy's 
dark matter halo. 

There is considerable disagreement in the existing lit- 
erature about Mvir,MV\^- Direct measures have typically 
focused on tracers of the inner halo, with recent results 
ranging from ~ [0.8 - 1.0] x 10^^ (Battaglia et al. 
2005; Smith et al. 2007; Xue et al. 2008) to 1.6 x lO^^ 
(Gnedin et al. 2010), while Deason ct al. (2012) argued 
that the mass of the halo within 150 kpc likely falls in the 
range of [0.5 — 1.0] x 10^^ Mq. Published measurements 
using the radial velocities of satellite galaxies give an 
even wider range, from ^ [0.8 — 2.5] x 10^^ Mq, depend- 
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ing on the sateUites included and assumptions about the 
orbits of those satellites (Kochanek 1996; Watkins ct al. 
2010). Boylan-Kolchin et al. (2011a) and Busha et al. 
(2011) combined analyses of large A^-body simulations 
with the Kallivayalil et al. (2006b) measurement of the 
Large Magellanic Cloud's proper motion and obtained 
Mvir.MW ~ 2 X 10^2 and 1.2 x 10^^ Mq, respectively. Pre- 
vious determinations of Mvir,MW based on the Leo I tim- 
ing argument have ranged from [1 — 3] x 10^^ Mq (Zarit- 
sky et al. 1989; Kochanek 1996; Li & White 2008, Paper 
I), while the Milky Way-M31 timing argument is also 
more consistent with values of Mvir,MW that are some- 
what higher than those obtained via stellar tracers (van 
der Marel et al. 2012b). 

Indirect constraints on Afvir,MW can be obtained from 
combinations of galaxy-galaxy lensing and Tully-Fisher 
data. Button et al. (2010) and Reyes et al. (2011) 
parametrize this combination of data in terms of the typ- 
ical ratio of a galaxy's optical velocity Vo-pt (the circular 
velocity at 2.2 disk scale radii) to its virial velocity^ V2ooc- 
Reyes et al. find Vo^t = l-27V2ooc for galaxies with stel- 
lar masses equal to that of the Milky Way (M^^mw ~ 
6 X W^Mq] Flynn et al. 2006; McMillan 2011); includ- 
ing scatter in concentration at fixed mass, this becomes 

^opt < l-8V'200c (2cr). Using Kpt .MW = 240 km s 
(McMillan 2011; Schonrich 2012), we find that Reyes et 
al. predict a median value of Mvh-.MW ~ 2.5 x 10^^ M©, 
with Mvir.MW > 1-16 X 10^2 Mq (2cr). This is simi- 
lar to abundance matching's prediction of Mvir,MW ~ 
2.5 X Mq for this stellar mass (Guo et al. 2010; 
Moster et al. 2012; Reddick et al. 2012). 

In light of these results, it is fair to say that the mass 
of the Milky Way is known to no better than a factor of 
two, and possibly even worse; indeed, Klypin ct al. (2002) 
showed that a wide range of observational data can be fit 
both for Mvir = 10^^ and for 2 x lO^'^. This uncertainty is 
a fundamental limitation in understanding several press- 
ing questions related to galaxy formation and cosmol- 
ogy. For example, it is frequently assumed that galaxies 
have available to them (fib/Om) x M^i^ in baryons. If 
-^vir,MW = 5 X 10^^ M0, then essentially all of those 
baryons are accounted for by observed stars and gas in 
and around the Galaxy; if Mvir, mw = 2.5 x 10^^ Mq, then 
a large majority of the Milky Way's baryons are missing 
(Fukugita & Peebles 2004; Anderson & Brcgman 2010; 
Gupta et al. 2012; Fang et al. 2012). Another exam- 
ple is the predicted abundance of dark matter satellites, 
which scales nearly linearly with dark matter halo mass; 
the expected mass of a galaxy's most massive satellite 
scales in the same manner. As a result, our interpre- 
tation of the missing satellites problem (Klypin ct al. 
1999; Moore ct al. 1999) and related questions about 
the Galaxy's satellite population depends on Mvir mw 
(Boylan-Kolchin ct al. 2011b, 2012; Vcra-Ciro et al. 2012; 
Wang ct al. 2012). 

A measurement of the mass distribution of the Milky 
Way at large Galacto-centric distance {D > 100 kpc) 
would bring significantly more clarity to each of these 
issues. Here, we use the results of Paper I in conjunction 

These authors define the virial velocity to be the halo's circular 
velocity at the radius containing an average density of 200 times 
the critical density. 



with cosmological simulations of Milky Way-size dark 
matter halos to provide a constraint on Mvir,MW based 
on the proper motion of Leo 1 [D = 260.6 kpc). This 
paper is structured as follows. Section 2 describes the 
simulations we employ and provides a summary of the 
results of Paper I. We then compare the simulations to 
the radial velocity (Section 3.1) and space velocity (Sec- 
tion 3.2) of Leo I. This allows us to constrain the Milky 
Way's virial mass (Section 4). We discuss our results in 
the context of expected infall times of ACDM satellites 
in general, and Leo I in particular, in Section 5, and con- 
clude with a brief summary and discussion of our results 
in Section 6. 

2. SIMULATIONS AND DATA 

2.1. Simulations 

Our primary set of Milky Way-size halos are those of 
the Aquarius simulation (Springcl ct al. 2008a). The 
simulations, and the properties of the halos, have been 
extensively discussed elsewhere (e.g., Springel et al. 
2008a, b; Navarro et al. 2010), and we refer the reader 
to those papers for details related to the simulations. 
We use the six halos (denoted A-F) simulated at "level 
two" resolution of ~ 10"* Mq per particle with Plummer- 
equivalent gravitational softening length of 66 pc. The 
Aquarius suite was run assuming a WMAP-1 cosmology 
with parameters fi™ — 0.25, fiA — 0.75, cs = 0.9, Hq — 
73kms^^ Mpc^^, Us — 1.0. More recent measurements 
prefer slightly different values of these parameters, but 
these differences are not expected to affect the assem- 
bly and internal dynamics of Milky Way-mass halos at a 
substantive level and are therefore unimportant for our 
analysis. In Appendix A, we compare results from sim- 
ilarly well-resolved halos - VL-II and GHALO (Stadel 
et al. 2009) - using WMAP3 parameters to show this 
explicitly. 

The virial mass Mvir of a dark matter halo is not 
uniquely defined, but rather depends on the density 
threshold used to define a dark matter halo. Throughout 
this work, we will define Mvir to be the mass within a 
sphere, centered on the halo in question, containing an 
average density A times the critical density of the Uni- 
verse; we use A = Avir, the value derived from the spher- 
ical top-hat collapse model (Gunn & Gott 1972; Bryan 
& Norman 1998), which gives Avir ~ 94 at z = for 
the cosmology of the Aquarius simulations. Note that 
other common choices of A are 200 and 20QQ.m{z). The 
Aquarius halos have 0.95 < Mvir/lO^^ Mq < 2.2, which 
is a reasonable approximation of current constraints on 
the value of the Milky Way's virial mass (see Boylan- 
Kolchin et al. 2012 for a recent compilation). 

In the ACDM cosmogony, all satellite galaxies are ex- 
pected to initially form within their own dark matter 
halo, outside of their host's virial radius. The luminous 
satellite galaxies of the Milky Way are therefore expected 
to be represented by some subset of dark matter subha- 
los surviving to z = 0. The precise relationship between 
satellite galaxies and dark matter subhalos is presently 
unclear; however, all models of galaxy formation predict 
a good correspondence between the most luminous satel- 
lite galaxies and the most massive subhalos for some suit- 
able definition of mass, most often the maximum mass 
or circular velocity that the sublialo has ever had (e.g.. 
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Kravtsov et al. 2004; Conroy et al. 2006; Guo ct al. 2010; 
Yang et al. 2012; Moster et al. 2012; Reddick ct al. 2012; 
Bchroozi et al. 2012). We therefore record the redshift 
at which a subhalo's bound mass is maximized, which we 
hereafter refer to as the infall redshift* ^mfaib as well as 
the maximum circular velocity at that time, VSntaii. 

Dark matter subhalos with maximum circular veloci- 
ties in excess of 5 kms~^ and having T4iax > 0.044 Vvir 
(corresponding to 5kms^^ for the lowest mass halo, Aq- 
B) are selected for further analysis. All subhalos within 
one Megaparsec of their host satisfying this Knax thresh- 
old are considered, with no additional velocity-based cri- 
teria. This ensures we do not introduce any bias against 
high velocity satellites in our analysis. All distances 
and velocities discussed below are computed in the halo- 
centric frame. 

2.2. Observational Data 

The observational data used in our analysis is detailed 
in Paper I. For the purposes of this work, the important 
quantities are: 

• The Galacto-centric distance to Leo I (I?looi): 

DlooI = 260.6 ± 13.3 kpc (1) 

• The Galacto-centric radial velocity of Leo I 

(Vr,Leol): 

Lcoi = 167.9 ± 2.8 km s"^ (2) 

• The Galacto-centric transverse (tangential) veloc- 
ity of Leo I (Vtxcoi): 

Vtxooi = 101.0 ±34.4 km s"^ (3) 

The measured three-dimensional Galacto-centric space 
velocity of Leo I is 

FlooI = 199.8+^^:^ [+^^;°] km s-i , (4) 

where the errors - computed from Monte Carlo sam- 
pling of the proper motion error space - represent 68.27% 
(95.45%) confidence intervals about the stated median. 

3. INTERPRETING THE MOTION OF LEO I 

3.1. The radial velocity of Leo I 

Figure 1 shows a phase space diagram for all subhalos 
in the Aquarius simulations. The radial velocity Vr of 
each subhalo is scaled by the virial velocity of its host, 
and the halo-centric radius of each subhalo is likewise 
scaled by the virial radius of the host. Inf ailing sub- 
halos have Vr < 0, while outgoing subhalos have posi- 
tive radial velocities. The black contours contain 68.3%, 
95.5%, and 99.7% of subhalos^. Additionally, we have 
placed square symbols with error bars^" on the plot to 
denote the observed position of Leo I in this phase space 

* This definition of infall time is similar to the time a subhalo 
last crossed the physical boundary of its host but does not rely on 
any specific definition of the host's "edge". 

^ The contour value at a specific radius for a given percentage 
value are computed by finding the location of that value in the 
cumulative distribution function of Vr at that radius. 

Unless otherwise noted, all plotted error bars are 1 a. 
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Fig. 1. — Radial phase space diagram of Aquarius subhalos, 
scaled to virial units. The inner, middle, and outer contours con- 
tain 68, 95.4, and 99.7% of subhalos from the Aquarius simulations. 
Colored squares with error bars show the location in phase space of 
Leo I, with 1 a uncertainties, for a variety of assumed virial masses 
for the Milky Way. Note that the absolute error on the radial 
velocity is smaller than the size of the symbols. The observed val- 
ues of (r, Vr) for Leo I are consistent with fewer than 0.3% of the 
Aquarius subhalos if the MW virial mass is less than 1.1 X lO"'^'^ A-Iq . 

for four representative values of the Milky Way's virial 
mass: Mvir,Mw/10^^ = 0.7 (black), 1.0 (gray), 1.5 
(magenta), and 2.0 (cyan). 

It is obvious why the observed radial velocity of Leo 
I has demanded a high virial mass for the Milky Way: 
subhalos in both the 7 x 10^^ and 10^^ Mq halos have 
radial velocities as high as Leo I less than 0.3% of the 
time. In more massive halos, Leo I's high radial velocity 
becomes more likely; however, it is still higher than 70% 
of subhalos at Leo I's Galacto-centric distance in a halo 
with Afvir = 2 X 10^^ Mq. With only a measurement 
of Leo I's radial velocity, the preferred mass of the MW 
is significantly in excess of 10^^ Mq. Constraining Mvir 
more tightly requires knowledge of the transverse veloc- 
ity of Leo I, however, as it requires understanding how 
likely it is for a subhalo to be on an unbound orbit about 
its host, and whether the observed tangential velocity of 
Leo I is cosmologically plausible, given its radial velocity. 
Furthermore, as we demonstrate in Section 5, satellites 
that have been recently accreted occupy a wide range 
of radial velocities but are confined to a much narrower 
range of total velocities. In combination with reasonable 
estimates of infall times (based on, e.g., star formation 
histories), 3D velocity information can therefore be much 
more constraining than radial velocity information alone. 

3.2. The space velocity of Leo I 

The analysis of section 3.1 does not require any knowl- 
edge of the transverse velocity of Leo I. In this section, 
we incorporate the measurement of Vt^coi from Paper I 
to compare the full space velocity of Leo I to simulated 
ACDM satellites. In Figure 2, we again plot the phase 
space diagram of Aquarius subhalos, this time using to- 
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Fig. 2. — The phase space diagram for Aquarius subhalos us- 
ing their space velocity; velocities and positions are scaled by 
virial values. At each radius, the velocity distribution is roughly 
Maxwellian rather than Gaussian; the distributions are maximized 
near Vtot ~ K/ir for ^ wide range of radii. The colored squares show 
Leo I's location in this phase space, using the measurement of Vt 
from Paper I, for representative values of M^ir.MW (the colors arc 
the same as in Fig. 1). Colored lines show surfaces of constant en- 
ergy in an NFW potential, while the dashed (dotted) curve shows 
the escape velocity for an NFW potential with a concentration of 
8 (16). For virial masses less than 1.2 X 10^^ A^0i Leo I is less 
bound than virtually all Aquarius subhalos. 

tal velocity V rather than radial velocity on the vertical 
axis. The direction of the radial velocity of each satellite 
is indicated by the sign of the total velocity: negative 
(positive) for inward (outward) radial velocities. 

We place Leo I on the plot for the same four choices 
of the Milky Way's virial mass used in Figure 1: 
Mvir,Mw/lO^^M0 = 0.7 (black), 1.0 (gray), 1.5 (ma- 
genta), and 2.0 (cyan). Lines of constant energy for an 
object with Leo Fs Galacto-centric distance and observed 
velocity for each of these halo masses are also plotted us- 
ing the same color scheme, assuming a dark matter halo 
following a Navarro, Frenk, & White (1997, hereafter 
NFW) profile with a concentration of c = 12 (Maccio 
et al. 2007). Also plotted in Figure 2 are curves corre- 
sponding to the escape velocity Vase as a function of ra- 
dius for an NFW profile with a concentration parameter 
of 8 (black dashed curve) and 16 (black dotted curve). In 
this representation, these curves are independent of the 
assumed halo mass. The values used (c = 8, 16) approx- 
imately span the expected 1 a scatter in concentration at 
fixed mass for Milky Way-size halos (Bullock et al. 2001; 
Maccio ct al. 2007; Boylan-Kolchin et al. 2010; Klypin 
et al. 2011). 

A remarkable feature of Fig. 2 is that vanishingly few 
subhalos - to be precise, only one in our entire Aquar- 
ius sample, or 0.01% - have orbits that are unbound 
(y > Vqsc)- Furthermore, this single unbound subhalo 
is part of a massive (Mvir ~ 1.3 x 10^^ Mq) infalling 
group^"'^ at r « 1.4rvir in the Aquarius C halo, and is 

This group can be seen as a deformation of the contours at 
Vr ~ — 1.5Vvir in Figure 1 



therefore a somewhat special case. Similar results are 
found for VL-II (0.06 % unbound) and GHALO (no un- 
bound subhalos) ; see Appendix A for a version of Fig 2 
that includes VL-II and GHALO data. In the ELVIS 
suite of simulated MW-size halos (Garrison-Kimmel et 
al., in preparation), which adopts a WMAP 7 cosmology, 
eight of ten halos have no unbound subhalos. Approx- 
imately 1% of subhalos in the other two halos, both of 
which have had a recent merger, are unbound. Dcason 
et al. (2011) and Di Cintio et al. (2012) find that approx- 
imately 3% of subhalos within 1 — l.Srvir are unbound; 
the origin of this difference is presently unclear. 

Based on our analysis of the eight highest resolution 
ACDM N-body simulations of Milky Way-sized dark mat- 
ter halos performed to date, it is very unlikely that Leo 
I is on an unbound orbit. This point is consistent with 
both theoretical models of halo formation and numerical 
simulations of structure formation, as discussed in the 
Introduction^^. Although unbound orbits are quite rare 
in the simulations adopted here, they are fully incorpo- 
rated in our analysis. Likewise, our analysis already in- 
cludes any subhalos with energies that have been boosted 
by 3-body interactions. That Leo I is almost certainly 
bound imposes a weak constraint on the MW mass in 
the context of ACDM, however: this requires only that 
-^^vir,MW ?i 0.7 X 10^^ M0 (ignoring proper motion er- 
rors).^'' In the next section, we combine the proper mo- 
tion results of Paper I with the Aquarius subhalo data 
introduced above to derive more stringent lower limits 
on the virial mass of the Milky Way. 

4. CONSTRAINING THE MILKY WAY'S VIRIAL MASS 

While the Aquarius halos provide us with a large sam- 
ple of subhalos, the host halos themselves only give us 
six different virial masses. The results of Figures 2 and 7 
indicate that halos in the mass range of interest for the 
Milky Way are very close to self-similar in terms of the 
kinematics of their subhalo populations, however, when 
each halo is scaled to virial quantities. We can therefore 
use our Aquarius sample to interpret Leo I's motion in 
a variety of halo masses. In practice, this simply means 
that to change from a halo with original virial quantities 
(M, R, V) to a halo with virial quantities (A/', R' , V), 
we multiply all radii by R' /R and all velocities by V /V. 

Given a subhalo's current position and velocity, we 
compute the velocity it would have at the Galactocentric 
distance of Leo I based on its binding energy. The or- 
bital energy of a subhalo is calculated using a spherically 
symmetric NFW profile with c = 12. We have compared 
this calculation to using the full gravitational potential 
computed from the parent A^-body simulation and find 
very good agreement, with a small level (< 10%) of sym- 
metric scatter at a given radius due to non-sphericity of 
the gravitational potential (see, e.g, Hayashi et al. 2007). 

The same logic can be applied to fast-moving satellites of M31, 
such as And XII and XIV (McConnachic 2012), meaning they are 
also likely to be bound satellites. 

Note that this highlights the importance of the mass distribu- 
tion exterior to Leo I's Galacto-centric distance: an NFW profile 
with Mvir = 7 X 10^^ and c = 8 places Leo I on a parabolic orbit 
(Vcsc = Vlooi)i whereas the same mass distribution, truncated at 
i^Leoli has an escape velocity that is 20% lower, ~ 160kms~^. 
The difference between extended and point-mass distributions is 
even more pronounced for higher values of Mvir- 
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The probability that a random subhalo i has a binding 
energy £ less than that of Leo I (i.e., that the subhalo 
is less bound than Leo I) - p{£i < Slcoi) ~ is then equal 
to the probability of subhalo i having a space velocity 
greater than that of Leo L We incorporate the Monte 
Carlo samplings of the proper motion error space de- 
scribed in Paper I when computing pi — p{£i < £looi)- 
For an ensemble of N subhalos, the probability that one 
subhalo chosen at random is less bound than Leo I is 

P{£ <£lcoi) = ^Y.P^- (5) 

i 

We consider only satellites with positive radial velocities 
(subhalos moving away from the halo center, as is the 
case for Leo I) in our analysis in order to make a fair 
comparison to the dynamics of Leo L 

The left panel of Figure 3 shows probability distri- 
butions for finding subhalos on orbits at least as ener- 
getic as that of Leo I. The dotted curve shows the re- 
sult when considering any individual galaxy. Even at 
Mvir = 2 X 10^^ M0, a randomly selected subhalo only 
has a 20% chance of being less bound than Leo L Of 
course, part of the motivation for obtaining measure- 
ments of Leo I's proper motion was its high radial ve- 
locity, i.e., Leo I is not a randomly chosen satellite. The 
probability of finding a high velocity satellite will obvi- 
ously increase as the sample size of satellites increases. 

Ideally, we would use proper motion data for a sta- 
tistically representative set of satellite galaxies. Unfor- 
tunately, only five Milky Way satellites have measured 
space velocities that are accurate at the 25% level^"*. This 
is still a data set without an obvious, homogeneous se- 
lection function, although we note it does consist of the 
five most luminous Milky Way satellites (excluding the 
Sagittarius). A blind search for Galactic satellites, ir- 
respective of luminosity, will preferentially find nearby 
satellites owing to luminosity bias. This bias does not 
affect searches within ^ 400 kpc of the Milky Way for 
satellites with Ly ^ 10^ L0, so long as there is not a 
population of very extended satellites above this lumi- 
nosity with surface brightnesses that fall below current 
detection limits. The 11 "classical" satellites of the Milky 
Way (e.g., Mateo 1998) therefore comprise a sample with 
a relatively well-known selection function. We can then 
ask how likely it is to find a satellite on an orbit as ener- 
getic as Leo I's when selecting from a sample of this size. 
Note that this should result in lower limits for the MW 
mass that are conservative (in the sense that different 
assumptions about the unknown transverse velocities of 
six classical satellites will only lead to an increased lower 
limit on Mvir,Mw)j as fewer than half of these 11 satellites 
have measured space velocities. 

The lack of strongly constraining proper motion data 
for 6 of the 11 classical satellites means that, in principle, 
there could be as many as 6 additional satellites on orbits 
as energetic as Leo I. We can estimate how likely this is 
on a satellite-by-satellite case by comparing to the sim- 
ulations by first selecting subhalos with similar Galacto- 
centric distance and observed radial velocity for a given 
satellite, then using the distribution of space velocities 

These satellites are the Large and Small Magellanic Clouds 
(Kallivayalil ot al. 2006b, a; Piatck ot al. 2008), Fornax (Piatck ct al. 
2007), Sculptor (Piatck ct al. 2006), and Leo I (Sohn ct al. 2012b). 



of these subhalos to compute the distribution of orbital 
energies. This allows us to estimate the probability that 
a given satellite is on an orbit more energetic than Leo I. 
This probability is very low - less than 1% - in each case. 
Among the five satellites with well-measured space veloc- 
ities, Leo I has the smallest binding energy. We note that 
using the measurements of Kallivayalil et al. (2006b), the 
LMC would have a binding energy that is comparable 
to, or even less than, that of Leo I. However, the revised 
value of the LMC's Galacto-centric velocity (Kallivayalil 
et al. 2012, submitted) places it on an orbit that is more 
bound than Leo I so long as Mvir.MW ^ 7 x 10^^ Mq. 
We therefore take the scenario in which Leo I is the least 
bound of the eleven classical satellites as our fiducial case. 

The solid black curve in the left panel of Figure 3 
shows the probability of precisely one subhalo from a 
randomly-selected sample of eleven subhalos having a 
binding energy equal to or less than that of Leo I. For 
comparison, we also show the probability of finding pre- 
cisely two galaxies from such a sample (solid gray curve). 
If we ask instead how often at least one subhalo has a 
binding energy equal to or less than that of Leo I, we 
obtain the dashed black curve in left panel Figure 3. 
This figure shows that Leo I's measured space veloc- 
ity argues against a low mass for the Milky Way (i.e., 
against Mvir,MW ^ 10^^ M©). It is vanishingly rare for 
a randomly-chosen subhalo to be as energetic as Leo I if 
A'^vir.MW < 10^^ Mq, and there is less than a 10% chance 
of finding at least 1 such subhalo when selecting a ran- 
dom sample of 11 from hosts with Mvir,MW < 10^^ Mq. 

While the left panel of Figure 3 shows the likelihood of 
finding subhalo(s) less bound than Leo I in a halo with 
a given Mvir, a potentially more interesting quantity is 
the posterior probability distribution of Afvir given Leo 
I's orbital energy, P{M^i,.\£). This can be easily calcu- 
lated using Bayes' theorem by combining the results of 
the left panel of Figure 3 with the prior probability dis- 
tribution P(A/vir) for the Milky Way mass. One natural 
choice of prior in ACDM is P(A/vir) ~ dn/dM^ir, i.e., 
the probability of a given halo mass is proportional to 
the halo mass function. Another logical ACDM-based 
prior is P(Mvir) = M^i,- dn / dM^i^ , which approximates 
a stellar mass weighting^^. We also consider a flat prior 
on Mvir, PiMvir) = constant. Since the mass function 
decreases monotonically with Mvir and is never shallower 
than dn/dMyir esc M~^^'^ (e.g., Boylan-Kolchin ct al. 
2009), the priors incorporating the mass function assign 
higher weights to low-mass halos than to those with high 
mass and will result in lower estimates of Afvir than the 
flat prior. 

The results of these calculations are presented in Ta- 
ble 1 and plotted in the right panel of Figure 3. We 
take as our fiducial estimate the case where Leo I is the 
most energetic galaxy out of the sample of the 11 classi- 
cal Milky Way satellites with a prior of Mvir rfn/dMvir; 
this results in a posterior probability distribution (PPD) 
for the Milky Way's virial mass with a median value of 
1.6 X 10^^ Mq and a symmetric 90% confidence interval of 
Mvir.MW G [1-0, 2.4] X 10^2 The PPD for this choice 

This weighting assumes that the probability distribution of 
host halo masses for a randomly selected star is proportional to 
the halos' stellar content, and is perhaps the most reasonable prior 
if the Sun can be considered a randomly-selected star. 
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Fig. 3. — Left: Probability of finding N galaxies less bound than Leo I as a function of the Milky Way virial mass Mvir when choosing 
among M galaxies with well-known energies. (TV, M) = (1, 1) is the probability of an individual galaxy being less bound than Leo I (dotted 
curve), while (1, 11) [solid black curve] and (2, 11) [solid gray curve] correspond to finding one or two such galaxies out of a sample of eleven. 
Dashed curves show the probability of having at least 1 (black) or 2 (gray) galaxies as energetic as Leo I. Right: Posterior probability 
distribution (PPD) of the Milky Way's virial mass for various assumptions about sample size and prior. Our fiducial case (shaded) is the 
estimate we consider most reliable, with Leo I being the most energetic of a sample of eleven satellite galaxies and adopting a prior for the 
virial mass based on the mass-weighted dark matter halo mass function. In this case, the median value of the PPD is Mvir = 1.6 X 10-*^^ 
with a 90% confidence interval of Mvir MW G [1-0, 2.4] X 10^^ Mq. If we instead assume that one additional MW satellite is more energetic 
than Leo I, we obtain a PPD given by the dashed curve, which is shifted to higher values of Mvir.MW. The dotted curve shows the PPD 
for the case where Leo I is the most energetic of 25 satellites (the approximate rmmber of known Milky Way satellite galaxies). Each of 
these PPDs assumes a prior that is proportional to the mass-weighted dark matter halo mass fimction, M^ir dn/ dM^ii ■ A robust result of 
our analysis is that Mvir MW > lO^'^ Mq at 95% confidence when considering the eleven classical satellites of the Milky Way, irrespective 
of the unknown tangential velocities of some of these satellites and of the choice of prior. 



TABLE 1 

Estimates of the mass of the Milky Way (in units of 10^^ M©). 





Vt =0 
(prior: dn/dM^ir) 


Measured Vt 
(prior: fiat) 


Measured Vt 
(prior: Mvir dn/dM^ir) 


Measured Vt 
(prior: dn/dM^ir) 


1 of 11 

2 of 11 

3 of 11 

4 of 11 


1.1, [0.7 - 1.6] 
1.3, [1.0 - 2.1] 
1.6, [1.1 - 2.7] 
1.9, [1.3 - 3.6] 


1.6, [1.1 - 2.5] 
2.1, [1.4 - 3.4] 
2.5, [1.7 - 4.8] 
3.1, [2.0 - 6.8] 


1.6, [1.0 - 2.4] 
2.0, [1.4 - 3.1] 
2.3, [1.6 - 4.1] 
2.8, [1.9 - 5.6] 


1.5, [1.0 - 2.2] 
1.9, [1.3 - 2.8] 
2.2, [1.6 - 3.6] 

2.6, [1.8 - 4.6] 



Note. — Column (1) indicates the number of satellites with binding energies as low as 
Leo I; in all cases, we assume the satellites are selected from a sample of eleven (the classical 
Milky Way satellites). Columns (2)-(5) show the median value and 90% confidence interval 
of the posterior probability distribution for the virial mass of the Milky Way, in units of 
10^2 Mq. Columns (2) assumes that Leo I has zero tangential velocity (3 cr away from the 
result of Paper I), whereas columns (3)-(5) adopt the measured tangential velocity and its 
error distribution but use different priors on P(Mvir). When using the measured tangential 
velocity, we find that Mvir,MW > 10^^ Mq at 95% confidence, a result that is independent 
of the choice of prior and the unknown tangential motions of some of the classical satellites. 
Our best estimate of Mvir MW assumes that Leo I is the most energetic classical satellite and 
uses a mass prior that is proportional to the mass-weighted dark matter halo mass function 
(row one, column four), giving a median value of Mvir,MW = 1-6 x lO^'^ Mq. 
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Fig. 4. — The probability that, when choosing from eleven sub- 
halos, the most energetic will have a total velocity less than 
vtot at DlboI- The probabilities are plotted for four values of 
A^vir,Mw/10^^ Mq: 0.7 (black), 1.0 (gray), 1.5 (magenta), and 
2.0 (cyan). The dashed black vertical line shows the measured Vr 
of Leo I (which is an absolute lower limit for Vlcoi)i while the solid 
gray vertical line and gray shaded region show Leo I's measured 
Vtot and the 68.3% confidence interval about this measurement. It 
is very unlikely to find the fastest-moving subhalo to have Vtot as 
high as Leo I unless Afvir.MW > 10^^ A^q. 

is plotted as the shaded gray region in the right panel 
of Figure 3. The median and 90% confidence interval of 
the PPD shift to larger values as the number of satellites 
as energetic as Leo I is allowed to increase: if one of the 
classical satellites is less bound than Leo I, then our fidu- 
cial estimates would change to Mvu-.mw = 1-9 x 10^^ Mq 
and Mvir,MW S [1.4,3.1] x lO^^ ^t 90% confidence. 
The choice of prior has a weaker effect, as can be seen 
from comparing values in a given row of Table 1. Even 
the effect of sampling from a larger number of galaxies is 
modest: even if we assume that all 25 known Milky Way 
dwarfs had measured proper motions and that Leo I was 
still the most energetic from this entire set, the median 
of the PPD would only decrease to 1.3 x 10^^ Mq with a 
90% confidence interval of Mvir,MW G [0.9, 1.7] x 10^^ 
for our fiducial prior. 

While median and upper limits of the PPD vary based 
on the choice of input parameters. Table 1 shows that 
the data constrain (at 95% confidence) the Milky Way 
mass to exceed 10^^ Mq, independent of prior or number 
of additional fast-moving satellites. This is the strongest 
and most robust constraint provided by the combination 
of numerical simulations and the Leo I proper motion 
measurement. Indeed, even if we assume that Leo I has 
no tangential motion, which is ruled out at 3 cr by the 
observations of Paper I, we still find that the median of 
the PPD for A/vir mw exceeds 10^^ Mq (see column 2 of 
Table 1). 

The results of Figure 3 can also be expressed in terms 
of the probability of observing an object with the velocity 
of Leo I for different values of Mvir,Mw; this is perhaps 
more intuitive, as it corresponds to the observed quanti- 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 



Fig. 5. — Similar to Figure 2, except each subhalo is colored 
according to its infall time (cosmic time, in Gyr). Subhalos that 
fell in long ago (yellow-orange) exhibit a wide range of orbital en- 
ergies, with the largest concentration of points at small radii and 
high binding energies. Recently accreted subhalos (dark red-black) 
occupy a much narrower range of orbital energies and are substan- 
tially less bound than the typical early-infalling subhalo. 

ties. In Figure 4, we plot the probability that the most 
energetic subhalo out of a sample of 11 subhalos has a ve- 
locity (at Dlcoi) less than V , as a function of V , for four 
different values of Mvir,MW- The observed radial velocity 
of Leo I is shown as a dashed black vertical line, while 
the observed total velocity of Leo I is shown as a solid 
gray vertical line; the shaded gray region shows the 68% 
confidence interval about Vtot- The 7 x 10^^ M^g Milky 
Way is ruled out at 95% confidence level by the radial 
velocity alone (see also Table 1); including the observed 
tangential velocity only strengthens this conclusion. In 
general, it is unlikely to find the most energetic satellite 
(from a sample of 11) moving at Leo Fs velocity unless 
1 < A/,ir,Mw/10i' Mq < 2. 

5. LEO I AND THE NATURE OF SATELLITE INFALL 

As discussed in the Introduction, we expect to find a 
relationship between a satellite's orbital energy and the 
time it was accreted onto its host. The existence of such 
a correlation could aid in a broader interpretation of Leo 
I's space velocity by constraining when Leo I fell into 
the Milky Way and whether or not it has made multiple 
pericentric passages. 

We explore the relationship between infall time and 
orbital energy in Figure 5, which shows a version of the 
total velocity phase space for Aquarius subhalos in which 
each subhalo is colored according to its infall time (mea- 
sured in cosmic time, with 13.7 Gyr being the present 
day). We clearly see an energy- infall time trend: while 
the early infalling objects (yellow-orange) cover a wide 
range of orbital energies, with the majority of subhalos 
being tightly bound and having apocenters of 0.75 Tvir or 
less, recently accreted subhalos (dark red-black) lie in a 
well-defined, narrow range of energies with correspond- 
ing apocenters of ss 2rvir. 

Recently accreted subhalos track a curve of constant 



The Mass of the Milky Way from Leo I's Motion 



9 



energy - corresponding to Wtot(''vir) ~ 1.15 Vvir - quite 
well, especially before they reach their first pericenters. 
This is fully consistent with the findings of previous 
ACDM-based analyses: Wetzel (2011) finds that the 
typical infall velocity of satellites for hosts of Afvir ~ 
10^2 M0 is 1.1 - 1.15 Kir (see also Benson 2005 and 
Khochfar & Burkert 2006; the excellent agreement of 
wtot(''vir) may be partially coincidental, however, as our 
definition of binding energy and virial velocity differ 
slightly from those of Wetzel). Figure 5 supports the 
orbit calculations of Paper I, which strongly favor the 
scenario in which Leo I fell into the MW within the last 
2 Gyr and has recently completed its first pericentric 
passage (see also Rocha et al. 2012). As noted in Paper 
I, this orbital history agrees very well with the observed 
star formation history of Leo I, which shows continuous 
star formation until < 1 Gyr ago, with enhancements at 

4.5 and 2 Gyr in the past, after which star formation 
ceased completely (T. Smecker-Hane et al., in prepara- 
tion). 

The confinement of first infall orbits to a narrow range 
of orbital energies shows how valuable measurements of 
transverse velocities can be for interpreting satellite dy- 
namics. We further emphasize this point in Figure 6 by 
focusing only on recently accreted subhalos and contrast- 
ing the resulting phase space using 3D velocities (left 
panel) with the phase space using only radial velocity 
information (right panel). The contrast is stark: while 
recently accreted subhalos occupy a well-defined and nar- 
row range of the 3D phase space, these satellites cover a 
wide range of radial velocities at every radius. Transverse 
velocities clearly add a great deal of information that is 
missing in the radial phase space. This also implies that 
subhalos are typically accreted with non-negligible an- 
gular momentum, and that radial orbit approximations 
result in substantial information loss. A further impor- 
tant point is that there are recently-accreted satellites 
with small radial velocities but high tangential veloci- 
ties, indicating that other Milky Way satellites may have 
binding energies similar to Leo I even if they have small 
measured radial velocities. 

The apocentric distance for galaxies with binding en- 
ergies consistent with that of Leo I in a halo of Afvir = 

1.6 X 10^^ Mq is ~ 650 kpc, comparable to the estimate 
of Peebles et al. (2011). The virial radius of such a halo 
is ^ 300 kpc, meaning this apocentric distance is simi- 
lar to the turn- around radius (Mamon ct al. 2004). The 
first apocenter of an orbit after turn-around is expected 
to be at ^ 90% of the turn-around radius in the SSIM 
(Bertschinger 1985; Ludlow ct al. 2009), so this large 
apocenter value is consistent with theoretical expecta- 
tions for a galaxy on its first infall. 

The results of Section 5 have interesting implications 
for satellites on first infall into their host galaxies: Fig- 
ure 5 shows that such satellites should populate a nar- 
row range of energies. For a measured distance and ra- 
dial velocity of a satellite, then, the only uncertainties 
are the halo mass and transverse velocity. Leo T pro- 
vides just such an example for the Milky Way: a dis- 
tant satellite {D = 407 kpc) with a low radial velocity 
{v,. = — 61kms^^ in the Galacto-centric frame), Leo T 
is the only known dwarf spheroidal / transition object 
near the Milky Way with substantial HI gas content. 



and is therefore almost certainly falling into the Milky 
Way for the first time. Our best-fitting virial mass of 
1.6 X 10^^ M0 leads to a prediction for the transverse 
velocity of Leo T of « 120 ± 20kms"\ i.e., Leo T's 
transverse velocity should be approximately two times 
larger than its measured radial velocity. Future mea- 
surements of Leo T's proper motion with HST (Program 
GO-12914; PI: T. Do) will verify or disprove this pre- 
diction. Proper motion measurements for other distant 
satellites such as Leo II {D = 235 kpc) and Canes Ve- 
natici I (D = 218 kpc) would also be of great interest. 

6. DISCUSSION 

While the radial velocity of Leo I has been the basis of a 
substantial body of work related to the mass of the Milky 
Way and the properties of its satellites, our measurement 
of the proper motion of Leo I (in Paper I) adds vital in- 
formation about Leo I's orbit. We have shown that it 
is a priori extremely unlikely for Leo I to be unbound 
to the Milky Way in ACDM, as vanishingly few subha- 
los have velocities exceeding the local escape velocity of 
their host in all of the A^-body simulations analyzed here. 
By itself, this is not strongly constraining. However, we 
have shown that the phase space structure of subhalos in 
numerical simulations is very regular. We have therefore 
coupled the proper motion measurement with the simu- 
lations and have presented a new method for combining 
these data sets to derive constraints on Mvii_MW- 

Our best estimate for Mvir.MW is 1.6 x 10"^^ Mq with a 
90% confidence interval of [1.6-2.4] x 10^^ While the 
precise central value and upper limit of this range depend 
somewhat on input assumptions, our strongest finding 
is A/vir.MW > 10^^ Mq at 95% confidence, independent 
of mass priors or kinematic data from additional satel- 
lites. Compared to previous determinations of Mvir.MW 
based on the Leo I timing argument, this best-estimate 
value and range are similar to Zaritsky et al. (1989) and 
lower than, but consistent with, Li & White (2008) and 
Sohn et al. (2012b). van der Marel et al. (2012b) have 
computed Afvir, mw based on the timing argument for 
the MW-M31 pair and find an average for Afvir, mw of 
1.63 X 10^^ Af0, in good agreement with our results. 

Our results disfavor "light" Milky Way models, 
wherein Afvir, mw ^ 10"^^ AfQ. An appeahng feature of 
such models is that they help mitigate issues (Klypin 
ct al. 1999; Moore ct al. 1999; Boylan-Kolchin et al. 
2011b, 2012) in reproducing observations of the Milky 
Way's satellites in the context of ACDM (Vera-Ciro et al. 
2012; Wang ct al. 2012; Starkenburg et al. 2012; Zolo- 
tov et al. 2012). Our work, however, shows that it is 
extremely difficult to reproduce the observed velocity 
of Leo I in these low mass models. While mass con- 
straints based on tracer stars in the Galactic halo have 
traditionally favored masses of ~ 10^^, we note that an 
NFW halo with Afvir = 1-6 x lO^^ ^^d c = 12 has 
Af(< 150 kpc) = 10^2 Af0, consistent with the BHB star 
constraint from Dcason et al. (2012). 

The results of this paper provide further motivation for 
measuring proper motions for all of the classical Galactic 

^® Lcpinc ct al. (2011) measured a transverse velocity of 265.2 it 
129.4kms~^ for Leo II; the central value is tantalizingly large, but 
is also consistent with zero at 2 cr. 
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Fig. 6. — Phase space plots in terms of total velocity (left panel) and radial velocity (right panel) for subhalos accreted within the last 4 
Gyr of cosmic time; the color scale is the same as in Figure 5. The curve of constant energy for an object with Leo I's observed velocity 
and Galacto-centric distance in a halo of Mvir mw = 1-5 x lO"*^^ jW"©, ''tot(j'vir) ~ 1-15 VVin is shown in magenta. The dashed and dotted 
curves in the left panel are the same as in Figure 2. Whereas recently accreted satellites occupy a narrow range of velocities at fixed radius 
(left), they span a wide range of radial velocities (right). 

reflect that the Galaxy's quiescent merger history over 
the past several Gyr is coming to an end. The Milky 
Way-M31 orbit, as inferred from recent measurement of 
M31's proper motion (Sohn et al. 2012a; van der Marel 
et al. 2G12a), ensures that the Milky Way's quiescent 
history has a maximum future duration of ^ 4 Gyr. 



satellites: Table 1 indicates that if just one more of these 
galaxies is found to have a binding energy as low as that 
of Leo I, then the 95% confidence value for the lower limit 
of Mvii-^MW would increase to « 1.3 x 10^^ AIq. Reduc- 
ing the proper motion errors for Leo II is an especially 
high priority because the current mean value of its space 
motion (Lcpine et al. 2011) is very large (266 km s^^ at a 
distance of 235 kpc, comparable to that of Leo I). More 
simulations are also needed to understand potential sys- 
tematics of the method presented here. In particular, 
simulations of Local Group - rather than Milky Way - 
analogs are needed if we are to truly capture the nature 
of the Milky Way satellite system. 

In addition to constraining the mass of the Milky Way's 
dark matter halo, the orbital analysis of Paper I and the 
AGDM simulations studied here both argue for a first 
infall scenario in which Leo 1 has only recently joined 
the Milky Way. Coupled with (1) convincing evidence 
that star formation persisted in Leo I until the past 0.5- 
1 Gyr (Smecker-Hane et al. 2009), and (2) the lack of HI 
detected in Leo I (Mri < 1.5 x 10^ Mq; Grcevich & Put- 
man 2009), this implies that gas in MW dwarf satellites 
can be either expelled or removed on time-scales shorter 
than one crossing time (see also Peebles et al. 2011). If 
the gas was expelled via internal processes, then further 
studies of Leo I may shed light on star formation feed- 
back. If the gas has been removed by ram pressure, Leo 
I may provide an interesting constraint on the density of 
hot gas in the Milky Way halo. 

A further implication of a recent Leo I infall, in con- 
junction with evidence that both Magellanic Clouds are 
also on their first infall (Besla et al. 2007; Boylan-Kolchin 
et al. 2011a), is that the Milky Way system is a dynamic 
one, with substantial late-time assembly. This active re- 
cent history at the dark matter halo level may initially 
seem odd in the context of the implied quiescent merger 
history for the Milky Way galaxy; however, it may simply 



ACKNOWLEDGMENTS 

We have enjoyed fruitful discussions with Frank van 
den Bosch, Alis Deason, Michael Kuhlen, Aaron Ludlow, 
and Tammy Smecker-Hane. The Aquarius Project is part 
of the program of the Virgo Consortium for cosmologi- 
cal simulations. We gratefully acknowledge the Aquar- 
ius, Via Lactea II, and GHALO collaborations for giving 
us access to their simulation data, and thank Michael 
Kuhlen for providing us with a subhalo catalog from the 
GHALO simulation and Mark Vogelsberger for providing 
the peculiar potential data for the Aquarius A halo. MB- 
K acknowledges support from the Southern California 
Center for Galaxy Evolution, a multi-campus research 
program funded by the University of California Office of 
Research. MB-K thanks the Kavli Institute for Theo- 
retical Physics and the organizers of the "First Galaxies 
and Faint Dwarfs" workshop for providing a stimulating 
environment during the development of this paper. GB 
acknowledges support from NASA through Hubble Fel- 
lowship grant HST-HF-51284.01-A. This work was sup- 
ported in part by the National Science Foundation un- 
der grants AST-1009973 and PHYll-25915. Support for 
this work was also provided by NASA through a grant 
for program GO-12270 from the Space Telescope Science 
Institute (STScI), which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., 
under NASA contract NAS5-26555. This research has 
made use of NASA's Astrophysics Data System. 



The Mass of the Milky Way from Leo I's Motion 



11 



REFERENCES 



Anderson, M. E., & Bregman, J. N. 2010, ApJ, 714, 320 
Ascasibar, Y., Hoffman, Y., & Gottlober, S. 2007, MNRAS, 376, 
393 

Battaglia, G. et al. 2005, MNRAS, 364, 433 
Behroozi, R S., Wechsler, R. H., & Conroy, C. 2012, 

arXiv:1207.6105 [astro-ph] 
Benson, A. J. 2005, MNRAS, 358, 551 
Bertschinger, E. 1985, ApJS, 58, 39 

Besla, G., Kallivayalil, N., Hernquist, L., Robertson, B., Cox, 
T. J., van der Marel, R. P., & Alcock, C. 2007, ApJ, 668, 949 

Boylan-Kolchin, M., Besla, G., & Hernquist, L. 2011a, MNRAS, 
414, 1560 

Boylan-Kolchin, M., Bullock, J. S., & Kaplinghat, M. 2011b, 

MNRAS, 415, L40 
— . 2012, MNRAS, 422, 1203 

Boylan-Kolchin, M., Springel, V., White, S. D. M., & Jenkins, A. 

2010, MNRAS, 406, 896 
Boylan-Kolchin, M., Springel, V., White, S. D. M., Jenkins, A., & 

Lcmson, G. 2009, MNRAS, 398, 1150 
Bryan, G. L., & Norman, M. L. 1998, ApJ, 495, 80 
Bullock, J. S., Kolatt, T. S., Sigad, Y., Somerville, R. S., 

Kravtsov, A. V., Klypin, A. A., Primack, J. R., & Dekel, A. 

2001, MNRAS, 321, 559 
Busha, M. T., MarshaU, P. J., Wechsler, R. H., Klypin, A., & 

Primack, J. 2011, ApJ, 743, 40 
Byrd, G., Valtonen, M., McCall, M., & Innanen, K. 1994, AJ, 

107, 2055 

Conroy, C, Wechsler, R. H., & Kravtsov, A. V. 2006, ApJ, 647, 
201 

Deason, A. J. et al. 2012, arXiv:1205.6203 [astro-ph] 
— . 2011, MNRAS, 415, 2607 

Dehnen, W., McLaughlin, D. E., & Sachania, J. 2006, MNRAS, 
369, 1688 

Di Cintio, A., Knebe, A., Libeskind, N. I., Hoffman, Y., Yepes, 

G., & Gottlober, S. 2012, MNRAS, 423, 1883 
Diemand, J., & Kuhlen, M. 2008, ApJ, 680, L25 
Diemand, J., Kuhlen, M., & Madau, P. 2007a, ApJ, 657, 262 
— . 2007b, ApJ, 667, 859 

Diemand, J., Kuhlen, M., Madau, P., Zemp, M., Moore, B., 

Potter, D., & Stadel, J. 2008, Nature, 454, 735 
Diemer, B., More, S., & Kravtsov, A. 2012, arXiv:1207.0816 

[astro-ph] 

Button, A. A., Conroy, C, van den Bosch, F. C, Prada, F., & 

More, S. 2010, MNRAS, 407, 2 
Fang, T., Bullock, J. S., & Boylan-Kolchin, M. 2012, ApJ, 

(submitted) 

Fich, M., & Trcmainc, S. 1991, ARA&A, 29, 409 

Fillmore, J. A., & Goldreich, P. 1984, ApJ, 281, 1 

Flynn, C, Holmbcrg, J., Portinari, L., Fuchs, B., & JahreiC, H. 

2006, MNRAS, 372, 1149 
Fukugita, M., & Peebles, P. J. E. 2004, ApJ, 616, 643 
Ghigna, S., Moore, B., Governato, F., Lake, G., Quinn, T., & 

Stadel, J. 1998, MNRAS, 300, 146 
Gill, S. P. D., Knebe, A.. & Gibson, B. K. 2005, MNRAS, 356, 

1327 

Gncdin, O. Y., Brown, VV. R., GcUer, M. J., & Kenyon, S. J. 

2010, ApJ, 720, L108 
Grcevich, J., & Putman, M. E. 2009, ApJ, 696, 385 
Gunn, J. E., & Gott, J. R. L 1972, ApJ, 176, 1 
Guo, Q., White, S., Li, C, & Boylan-Kolchin, M. 2010, MNRAS, 

404, 1111 

Gupta, A., Mathur, S., Krongold, Y., Nicastro, F., & Galeazzi, 

M. 2012, ApJ, 756, L8 
Hayashi, E., Navarro, J. F., & Springel, V. 2007, MNRAS, 377, 50 
Kahn, F. D., & Woltjer, L. 1959, ApJ, 130, 705 
Kallivayahl, N., van der Marel, R. P., & Alcock, C. 2006a, ApJ, 

652, 1213 

Kallivayahl, N., van der Marel, R. P., Alcock, C., Axelrod, T., 
Cook, K. H., Drake, A. J., & Geha, M. 2006b, ApJ, 638, 772 
Khochfar, S., & Burkert, A. 2006, A&A, 445, 403 
Klypin, A., Kravtsov. A. V., Valenzuela, C, & Prada, F. 1999, 
ApJ, 522, 82 

Klypin, A., Zhao, H., & Somerville, R. S. 2002, ApJ, 573, 597 



Klypin, A. A., Trujillo-Gomez, S., & Primack, J. 2011, ApJ, 740, 
102 

Kochanek, C. S. 1996, ApJ, 457, 228 
Komatsu, E. et al. 2011, ApJS, 192, 18 

Kravtsov, A. V., Gnedin, O. Y., & Klypin, A. A. 2004, ApJ, 609, 
482 

Lepine, S., Koch, A., Rich, R. M., & Kuijken, K. 2011, ApJ, 741, 
100 

Li, Y.-S., & White, S. D. M. 2008, MNRAS, 384, 1459 
Ludlow, A. D., Navarro, J. F., Springel, V., Jenkins, A., Frenk, 

C. S., & Helmi, A. 2009, ApJ, 692, 931 
Maccio, A. V., Dutton, A. A., van den Bosch, F. C, Moore, B., 

Potter, D., & Stadel, J. 2007, MNRAS, 378, 55 
Mamon, G. A., Sanchis, T., Salvador-Sole, E., & Solanes, J. M. 

2004, A&A, 414, 445 
Mateo, M., Olszewski, E. W., & Walker, M. G. 2008, ApJ, 675, 

201 

Mateo, M. L. 1998, ARA&A, 36, 435 

McConnachic, A. W. 2012. AJ, 144, 4 

McMillan, P. J. 2011, MNRAS, 414, 2446 

Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., 

Stadel, J., & Tozzi, P. 1999, ApJ, 524, L19 
Moster, B. P., Naab, T., & White, S. D. M. 2012, 

arXiv: 1205.5807 [astro-ph] 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 

493 

Navarro, J. F. et al. 2010, MNRAS, 402, 21 
Peebles, P. J. E., Tully, R. B., & Shaya, E. J. 2011, 

arXiv:1105.5596 [astro-ph] 
Piatek, S., Pryor, C, Bristow, P., Olszewski, E. W., Harris, H. C, 

Mateo, M., Minniti, D., & Tinney, C. G. 2006, AJ, 131, 1445 
— . 2007, AJ. 133, 818 

Piatek, S., Pryor, C, & Olszewski, E. W. 2008, AJ, 135, 1024 
Reddick, R. M., Wechsler, R. H., Tinker, J. L., & Behroozi, P. S. 

2012, arXiv:1207.2160 [astro-ph] 
Reyes, R., Mandelbaum, R., Gunn, J. E., Pizagno, J., & Lackner, 

C. N. 2011, MNRAS, 417, 2347 
Rocha, M., Peter, A. H. G., & Bullock, J, 2012, MNRAS, 425, 231 
Sales, L. V., Navarro, J. F,, Abadi, M. G., & Steinmetz, M. 2007, 

MNRAS, 379, 1475 
Schonrich, R. 2012, arXiv:1207.3079 [astro-ph] 
Smecker-Hane, T. A., Marsteller, B., Cole, A., Bullock, J., & 

Gallagher, J. S. 2009, in Bulletin of the American Astronomical 

Society, Vol. 41, #419.11 
Smith, M. C. et al. 2007, MNRAS, 379, 755 
Sohn, S. T., Anderson, J., & van der Marel, R. P. 2012a, ApJ, 

753, 7 

Sohn, S. T., Besla, G., van der Marel, R. P., Boylan-Kolchin, M., 

Bullock, J. S., & Majewski, S. R. 2012b, ApJ, (submitted) 
Sohn, S. T. et al. 2007, ApJ, 663, 960 
Springel, V. et al. 2008a, MNRAS, 391, 1685 
— . 2008b, Nature, 456, 73 

Stadel, J., Potter, D., Moore, B., Diemand, J., Madau, P., Zemp, 

M., Kuhlen, M., & Quilis, V. 2009. MNRAS, 398, L21 
Starkcnburg, E. et al. 2012, arXiv: 1206.0020 [astro-ph] 
van der Marel, R. P., Besla, G., Cox, T. J., Sohn, S. T., & 

Anderson, J. 2012a, ApJ, 753, 9 
van der Marel, R. P., Fardal, M., Besla, G., Beaton, R. L., Sohn, 

S. T., Anderson, J., Brown, T., & Guhathakurta, P. 2012b, 

ApJ, 753, 8 

Vera-Ciro, C. A., Helmi, A., Starkcnburg, E., & Breddels, M. A. 

2012, arXiv: 1202.6061 [astro-ph] 
Wang, H., Mo, H. J., & Jing, Y. P. 2009, MNRAS, 396, 2249 
Wang, J., Frenk, C. S., Navarro, J. F., Gao, L., & Sawala, T. 

2012, MNRAS, 424, 2715 
Watkins, L. L., Evans, N. W., & An, J. H. 2010, MNRAS, 406, 

264 

Wetzel, A. R. 2011, MNRAS, 412, 49 
Xue, X. X. et al. 2008, ApJ, 684, 1143 

Yang, X., Mo, H. J., van den Bosch, F. C, Zhang, Y., & Han, J. 

2012, ApJ, 752, 41 
Zaritsky, D., Olszewski, E. W., Schommer, R. A., Peterson, 

R. C, & Aaronson, M. 1989, ApJ, 345, 759 
Zolotov, A. et al. 2012, arXiv:1207.0007 [astro-ph] 



12 



M. Boylan-Kolchin et al. 




Fig. 7. — Same as Fig. 2, but also including simulations with WMAP3 parameters (VL-II: orange circles; GHALO: blue triangles). 
Subhalos from these simulations populate the phase space diagram in an identical manner to the Aquarius subhalos; our Aquarius-based 
results should therefore be robust to modest changes in cosmological parameters. Note that VL-II has a mass of Mvir = 1-7 x 10^'^ Mq, 
which is near the upper end of the Aquarius halo masses, while GHALO has A/vir = 1.08 X 10^^ A^q, which is near the lower end. The 
virial scalings used throughout this work appear to be appropriate. 



APPENDIX 

COSMOLOGICAL PARAMETER (IN)DEPENDENCE 

A possible concern in interpreting the orbit of Leo I through the use of the Aquarius simulations is that the 
cosmological parameters adopted for Aquarius differ slightly from the currently favored values. Specifically, WMAP7 
results indicate that = 0.816 ± 0.024, fl^ = 0.274 ± 0.011, and = 0.968 ± 0.012 (Komatsu et al. 20il), placing 
the Aquarius parameters 3-10% off of the most recently measured values. To investigate the effects of variations in 
cosmological parameters. Figure 7 duplicates Figure 2 but also includes data for subhalos from two simulations using 
WMAP3 parameters {^m = 0.237, as = 0.742, = 0.951): VL-II (orange circles) and GHALO (blue triangles). 
Together, first and third-year WMAP parameters bracket parameters determined from the seven-year WMAP data 
release. The subhalos from WMAP3-based simulations populate phase space identically to those from the Aquarius 
simulations, indicating that small changes to cosmological parameters will have no effect on interpretations of Leo I's 
motion. 



